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We investigate the origin of the regular features observed in numerical studies of the interacting 
boson model with random interactions, in particular the dominance of L = ground states and the 
occurrence of vibrational and rotational band structures. It is shown that all of these properties 
can be interpreted and explained in terms of a Hartree-Bose mean-field analysis, in which different 
■ regions of the parameter space are associated with geometric shapes. The same conclusions hold for 

' the vibron model. 
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, O ' In empirical studies of medium and heavy even-even nuclei very regular features have been observed, such as the 
, tripartite classification of nuclear structure into seniority, anharmonic vibrator and rotor regions Q . In each of these 
three regimes, the energy systematics is extremely robust, and the transitions between different regions occur very 
rapidly, typically with the addition or removal of only one or two pairs of nucleons. Traditionally, this regular behavior 
has been interpreted as a consequence of particular nucleon-nucleon interactions, such as an attractive pairing force 
^ ' in semimagic nuclei and an attractive neutron-proton quadrupole-quadrupole interaction for deformed nuclei. 
On . It came as a surprise, therefore, that recent studies of even-even nuclei in the nuclear shell model (|-^ and in 
the interacting boson model (IBM) |^|-^ with random interactions (centered around zero, i.e. equally likely to be 
^—^ attractive or repulsive) displayed a high degree of order. Both models showed a marked statistical preference (> 60 %) 
I for ground states with L — 0, despite the random nature of the interactions. In addition, in the shell model evidence 
[ was found for the occurrence of pairing properties Q, and in the IBM for both vibrational and rotational band 
. structures 

' However, these results were obtained from numerical studies. It is thus important to gain a better understanding 
, as to why this happens. Can these properties be explained in a more intuitive way? Several attempts have been 
' made to understand these surprising phenomena, most of them for fermionic systems and especially with regards to 
O the dominance of L = ground states. Among others, we mention a study of the time-reversal invariance of random 
interactions , the connection with the width of the energy distributions |^J^ , the distribution of the lowest energy 
eigenvalues , the effect of higher-order interactions , the geometric chaoticity of the angular momentum coupling 
E> ■ of individual particles |^ , the connection with random polynomials , the correlation between wave functions and 
'k^ [ energies [ pT[ , the probability distribution of matrix elements p2[ , the relation with the diagonal matrix elements of 
. the Hamiltonian |l^ , and average energies and variances |0 . 

; In I it was shown that the overlap between the ground state wave function obtained in a single-j shell model 
calculation with random interactions with the seniority zero state is small, and that the distribution of these overlaps 
follows closely the predictions for chaotic wave functions. This suggests that for fermion systems, at least for like- 
nucleons in a single-j shell, the dominance of L = ground states may be due to the geometric chaoticity of randomly 
coupled individual spins, although the observed oscillations with shell size cannot be explained by such a simple 
scheme. On the other hand, in the IBM, besides the L = ground state dominance, a strong preponderance of 
vibrational and rotational values for energy ratios of yrast states was found, as well as a strong correlation with the 
corresponding vibrational and rotational values of the quadrupole transitions. This has been interpreted as evidence 
for the existence of vibrational and rotational structure [^|j6| . It is the purpose of this Rapid Communication to study 
the origin of the regular features observed in the IBM, when the parameters are chosen randomly. 

Collective excitations in nuclei are described in the IBM in terms of a system of interacting bosons [p^ . Its building 
blocks are a quadrupole boson with = 2+ and a scalar boson s'' with = 0+. The total number of bosons N 
is conserved by the IBM Hamiltonian. We consider the most general one- and two-body Hamiltonian 
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where the nine coefficients (eg, Cd, uq, U2, cq, C2, C4, uq, W2) are chosen independently from a Gaussian distribution of 
random numbers with zero mean and width a ||,|. 

The connection between the IBM, its potential energy surfaces, equilibrium configurations and geometric shapes, 
can be studied with mean- field Hartree-Bose techniques by means of coherent states 16,1^]. Since the one- and 



two-body IBM Hamiltonian of Eq. (g) does not give rise to a triaxial deformation [g^, the coherent state can be 
written as an axially symmetric condensate 

\N,a) = ( cos as^ + sin a dl ) |0) , (3) 



with — 7r/2 < a < Tr/2. The angle a is related to the deformation parameters in the intrinsic frame, /3 and 7 ]l5| , p^ . 
The potential energy surface is then given by the expectation value of the Hamiltonian in the coherent state 

En (a) — 04 sin^ a -f 03 sin'^ a cos a + 02 sin^ a + oq , (4) 

where the coefficients are linear combinations of the parameters of the Hamiltonian 
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ao = es + -Uq . (5) 

For random interactions, the trial wave fimction \N,a) and the corresponding energy surface En (a) provide infor- 
mation on the distribution of shapes that the model can acquire. The equilibrium configuration, or the intrinsic 
vibrational state, is characterized by the value of a = ao for which the energy surface attains its minimum value (this 
procedure is equivalent to solving the Hartree-Bose equation). For a given Hamiltonian, the value of uq depends on 
the coefficients a^, 03 and 02. Hence, the distribution of shapes for an ensemble of Hamiltonians depends on the joint 
probability distribution P{a4, a^, 02), which for the present case is given by a multivariate normal distribution. In 
practice, for each Hamiltonian the minimum of the energy surface E{a) is determined numerically. The equilibrium 
configurations can be divided into three different classes: an s-boson or spherical condensate {ao = 0), a deformed 
condensate with prolate or oblate symmetry (0 < ao < '"'/S or —tt/2 < ao < 0, respectively), and a d-boson condensate 
(ao-7r/2).^^ 

Each equilibrium configuration has its own characteristic angular momentum content. Even though we do not 
explicitly project the angular momentum states from the coherent state, the angular momentum of the ground state 
can be obtained from the rotational structure of the condensate in combination with the Thouless-Valatin formula 
for the corresponding moments of inertia. This procedure is described in detail in JlTf . The results are summarized 
in Table |. 

The s-boson condensate corresponds to a spherical shape. Whenever such a condensate occurs (in 39.4 % of the 
cases), the ground state has L — 0. 



• The deformed condensate corresponds to an axially symmetric deformed rotor. The ordering of the rotational 
energy levels L = 0, 2, . . . , 27V is determined by the sign of the moment of inertia 

Erot = ^^HL + 1) . (6) 

The deformed condensate occurs in 36.8 % of the cases. For X3 > the ground state has L = (23.7 %), while 
for X3 < the ground state has the maximum value of the angular momentum L = 2N (13.1 %). 
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• The d-boson condensate corresponds to a quadrupole oscillator with N quanta. Its rotational structure has a 
more complicated structure than the previous two cases. It is characterized by the labels t, and L. The 
boson seniority r is given by r = Stt-a + A = A^, — 2, . . . , 1 or for N odd or even, and the values of the 
angular momenta are L = A, A + 1, . . . , 2A — 2, 2A In this case, the rotational excitation energies depend on 
two moments of inertia 

Erot = ^T(r + 3) + ^L(L + 1) , (7) 

which are associated with the spontaneously broken three- and five-dimensional rotational symmetries of the 
d-boson condensate. The c?-boson condensate occurs in the remaining 23.8 % of the cases. The results in Table || 
can be understood qualitatively as follows. For > the ground state has r = for TV even or r = 1 for 
odd (~ 4 %), while for X5 < the ground state has the maximum value of the boson seniority r = iV (~ 19 
%). For T = and r = 1 there is a single angular momentum state with L = and L = 2, respectively. For the 
T — N multiplet, the angular momentum of the ground state depends on the sign of the moment of inertia I3. 
For X3 > the ground state has L = for iV = ik ox L = 2 for iV 7^ 3A; (9 %), while for X3 < the ground state 
has the maximum value of the angular momentum L = 2N (10 %). 

Table | shows that the spherical and deformed condensates contribute constant amounts of 39.4 % and 23.7 %, 
respectively, to the L = ground state percentage, whereas the contribution from the d-boson condensate depends 
on the number of bosons TV. The L = 2 ground states arise completely from the d-boson condensate solution. In 
Fig. we show the percentages of ground states with L — and L = 2 as a function of the total number of bosons 
N. A comparison of the results of the mean-field analysis (dashed lines) and the exact ones (solid lines) shows a good 
agreement. There is a dominance of ground states with L ^ for ^ 63-77 % of the cases. Both for L = and L — 2 
there are large oscillations with A^, which are entirely due to the contribution of the d-boson condensate. For N = 3k 
we see an enhancement for L = and a corresponding decrease for L = 2. In the mean-field analysis, the sum of the 
two, which corresponds to ~ 77 % of the cases, hardly depends on the number of bosons, in agreement with the exact 
results. For the remaining 23 % of the cases, the ground state has the maximum value of the angular momentum 
L = 2N. 

For the cases with an L = ground state, the probability distribution P{R) of the energy ratio 

E21 — Eq^ 

can be used to study the characteristics of the energy spectra. As mentioned above, in a numerical study of the IBM 
with random interactions Q it was found that P{R) exhibits two very pronounced peaks, right at the vibrational value 
of i? = 2 and at the rotational value oi R = 10/3. This provides a clear indication of the occurrence of vibrational 
and rotational structure, which was confirmed by a simultaneous study of the quadrupole transitions between the 
levels In Figs. || and y we show the contribution of each one of the equilibrium configurations to P{R) for A^ = 16 
and A^ — 15, respectively. For both cases, the spherical shape (solid line) contributes almost exclusively to the peak 
at i? = 2, and similarly the deformed shape (dashed line) to the peak at i? = 10/3, which once again confirms the 
vibrational and rotational character of these maxima. For A^ = 16 the contribution of the d-boson condensate (dotted 
line) is small, whereas for A^ = 15 it gives a contribution for small values of R, which corresponds to a level sequence 
L = 0, 4, 2. 

The above analysis is not only valid for the IBM, but can also be applied to other many-body systems of randomly 
interacting bosons. As an example, we mention the vibron model which was introduced to describe the relative 
motion in two-body problems, e.g. diatomic molecules |^|, nuclear clusters [ pO[ and mesons [ pl| . The vibron model 
has many of the same qualitative features as the IBM, namely vibrational and rotational spectra, but has a much 
simpler mathematical structure. Its building blocks are a dipole boson p'' with = 1~ and a scalar boson with 
— 0+. The total number of bosons A^ is conserved by the Hamiltonian. In this case, the coherent state is expressed 
as a condensate of a superposition of a scalar and a dipole boson j2^ 

\N,a) — I cos as'' -I- sin a pg I |0) , (9) 



with < a < Tr/2. The potential energy surface of a one- and two-body vibron Hamiltonian is a quadratic function 
of sin^ a 

Ejsi{a) = 04 sin'* a -|- 02 sin^ a -|- oq . (10) 
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In this case, the value of that characterizes the equihbrium configuration only depends on the two coefficients 04 
and a2. The distribution of shapes for an ensemble of Hamiltonians is determined by the joint probability distribution 
P{a4,a2). For the vibron model, the structure of the energy surface is simpler than for the IBM. This makes it 
possible to obtain analytic results for the probability of the occurrence of a given equilibrium configuration. As for 
the IBM, the solutions can be divided into three different classes. 

• The s-boson or spherical condensate (ao — 0) occurs with probability 

1 f f27\ 
Pi = — TT + 2 arctan W — = 0.396. (11) 





This solution corresponds to a spherical shape which has L = 0. 

• The deformed condensate (0 < ao < 7r/2) has probability 

1 /64 

P2 = — arctanW— = 0.216. (12) 
27r V 3 

The ordering of the rotational energy levels i = 0, 1, . . . , is determined by the sign of the moment of inertia 

Erot = • (13) 

For X3 > the ground state has L ^ (13.8 %), while for X3 < the ground state has the maximum value of 
the angular momentum L = N (7.8 %). 

• The p-boson condensate (ao = '^/'^) ha-s probability 

P3 = I-P1-P2 = 0.388 . (14) 



This solutions corresponds to a three-dimensional harmonic oscillator with N quanta. Its angular momentum 
content is given by L = A^, — 2, . . . , 1 or for A^ odd or even, respectively. For X3 > the ground state has 
L = for A^ even or L = 1 for A^ odd (17.9 %), while for X3 < the ground state has the maximum value of 
the angular momentum L = N (20.9 %). 

Table || shows that the spherical and deformed condensates contribute a constant amount of 53.4 % to the L — 
ground state percentage. For the p-boson condensate, the angular momentum of the ground state depends on the 
sign of the moment of inertia and the number of bosons A^. For even values of A^ , the ground state either has L = 
01 L = N, whereas for odd values it has either L = 1 or L = A^. As a result, the vibron model shows a dominance of 
L = ground states which oscillates between ^ 53 % for odd values of A^ and ~ 71 % for even values. The L = 1 
ground states arise completely from the p-boson condensate solution. In Fig. ^ we show the percentages of ground 
states with L — and L = 1 as a function of the total number of bosons A^. A comparison of the results of the 
mean-field analysis (dashed lines) and the exact ones (solid lines) shows an excellent agreement. The oscillations in 
the percentages of L = and L — 1 ground states are entirely due to the contribution of the p-boson condensate, 
whereas their sum is a constant which does not depend on A^. 

In this paper, we have studied the properties of low-lying states in the IBM with random interactions. We addressed 
the origin of the regular features, that were obtained before in numerical studies of the IBM, in particular the 
dominance of L = ground states and the occurrence of vibrational and rotational band structures. We have shown 
that a mean-field analysis can explain these features and provide a more intuitive understanding of their origin. The 
use of mean-field techniques and coherent states circumvents the use of coefficients of fractional parentage, bypasses 
the diagonalization of thousands of matrices, and makes it possible to associate different regions of the parameter 
space with particular intrinsic vibrational states, which in turn correspond to definite geometric shapes. For the IBM, 
there are three different equilibrium configurations or geometric shapes: a spherical shape (~ 39 %), a deformed 
shape (~ 36 %) and a condensate of quadrupole bosons (^ 25 %). Since the spherical shape only has L = 0, and the 
deformed shape has this value of L in about two thirds of the cases, these two solutions account for ~ 63 % of L = 
ground states. The oscillations observed for the L = ground state percentage can be ascribed completely to the 
contribution of the c?-boson condensate. In addition, we found a one-to-one correspondence between the occurrence 
of the spherical and deformed equilibrium configurations and the peaks in the probability distribution for the energy 
ratio R — 2 and R = 10/3, respectively. Similar conclusions hold for the vibron model. 
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The study of interacting boson models with random interactions indicate that there is a significantly larger class 
of Hamiltonians that leads to regular, ordered behavior at the low excitation energies than was commonly assuemd. 
The fact that these properties are shared by both the IBM and the vibron model, seems to exclude an explanantion 
solely in terms of the angular momentum algebra, the connectivity of the model space, or the many-body dynamics 
of the model, as has been suggested before. The present analysis points to a more general phenomenon that does not 
depend so much on the details of the angular momentum coupling, but rather on the occurrence of definite, robust 
geometric phases such as spherical and deformed shapes. 

Since the collective subspace involves a drastic truncation of the shell model space as well as a bosonic approximation, 
our conclusions cannot be applied directly to the fermion space. We are currently exploring a nucleon-pair truncation 
scheme in order to search for answers in this more complex situation p3[ . These results may help us understand the 
appearance of robust properties in many-body quantum systems with random interactions, be it in nuclear, atomic, 
molecular or mesoscopic systems like quandom dots [Q. 
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TABLE I. Percentages of ground states with L = 0, 2 and 2N , obtained in a mean-field analysis of the interacting boson 
model. 
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TABLE IL Percentages of ground states with L 


= 0, 1 and N, obtained in 


a mean-field analysis 


of the vibron model. 
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FIG. 1. Percentages of ground states with L = and L = 2 in the IBM with random one- and two-body interactions 
calculated exactly for 10,000 runs (solid lines) and in mean- field approximation (dashed lines). 
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FIG. 2. Probability distribution P{R) of the energy ratio R obtained for N = 16 and 10,000 runs for the spherical (solid 
line), deformed (dashed line) and d-boson condensate (dotted line) equilibrium configurations, respectively. 
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